Macroscopic magnetization jumps due to independent magnons 
in frustrated quantum spin lattices 
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For a class of frustrated spin lattices including the kagome lattice we construct exact eigenstates 
consisting of several independent, localized one-magnon states and argue that they are ground states 
for high magnetic fields. If the maximal number of local magnons scales with the number of spins 
in the system, which is the case for the kagome lattice, the effect persists in the thermodynamic 
limit and gives rise to a macroscopic jump in the zero-temperature magnetization curve just below 
the saturation field. The effect decreases with increasing spin quantum number and vanishes in the 
classical limit. Thus it is a true macroscopic quantum effect. 
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In frustrated quantum spin lattices the competition of 
quantum and frustration effects promises rich physics. 
A reliable description of such systems often constitutes a 
challenge for theory. A famous example is the kagome lat- 
tice antiferromagnet. In spite of extensive studies during 
the last decade its ground state properties are not fully 
understood yet. Classically it has infinite continuous de- 
generacies. In the quantum case (s=l/2), the system is 
likely to be a spin liquid with a gap for magnetic excita- 
tions and a huge number of singlet states below the first 
triplet state (see jl[ || [| and references therein). 

In this Letter we will focus on the zero-temperature 
magnetic behavior of highly frustrated lattices, in par- 
ticular for high magnetic fields. One aspect is given by 
the observation of nontrivial magnetic plateaus in frus- 
trated two dimensional (2D) quantum antiferromagnets 
like SrCu2(B03) j|, 0], which has stimulated theoreti- 
cal interest (see e.g. ||). Also the kagome lattice has a 
plateau at one third (m = 1/3) of the saturation magne- 
tization ^ . Since this plateau can be found also in the 
Ising model and in the classical Heisenberg model with 
additional thermal fluctuations [H it can be considered 
to be of classical origin. However, the structure of the 
ground state in the classical model is highly non-trivial 
at m — 1/3 pJ] and has not been clarified yet for the 
quantum model. 

Another aspect is given by unusual jumps seen in mag- 
netization curves. Such jumps can arise for different 
reasons. One possibility is a first-order transition be- 
tween different ground states like the spin flop transition 
in classical magnets or in strongly anisotropic quantum 
chains jO]. Here we discuss another possibility, namely 
a macroscopically large degeneracy in the exact ground 
states of the full quantum system for a certain value of 
the applied field. We argue that this is a general phe- 
nomenon in highly frustrated systems. This is remark- 
able in so far as one can exactly write down ground states 
at a finite density of magnons in a strongly correlated 
system which is neither integrable, nor has any apparent 



non-trivial conservation laws. Such jumps represent a 
genuine macroscopic quantum effect which is also of pos- 
sible experimental relevance since it occurs in many well- 
known models like the kagome lattice. This jump occurs 
just below saturation and should be observable in mag- 
netization experiments on the corresponding compounds 
if the coupling constants are small enough to make the 
saturating field accessible. 

We consider N quantum spins of "length" s described 
by the Hamiltonian 



H = ^2 J ij 

m 
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where the sum runs over neighboring sites (ij), S z is 
the z-component of the total spin S z — Y^i^ii h the 
magnetic field, A the AAZ-anisotropy and are the 
exchange constants. 

If the magnetic field h is sufficient large, the ground 
state of (Q) becomes the magnon vacuum state |0) = 
| ttt . . . ttt) where all spins assume their maximal Sf- 
quantum number. The lowest excitations for the case of 
a high magnetic field are one-magnon states |1). They 
are represented by states where the 5" z -quantum number 
is lowered by one and can be written as 



1 N 

i) = -53oisrio), 
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where c is chosen such that (1|1) = 1. For a lattice with 
n s spins per unit cell, there are n s possible magnon bands 
Wi{k) with i = 0, . . . , (n s — 1). For certain combinations 
of Jij of the Hamiltonian ([!]) the lowest magnon disper- 
sion Wo(k) becomes flat (i.e. independent of k) in some 
directions or in the whole /c-space. 

If the one-magnon dispersion is independent of one of 
the components of k, one can use Fourier transformation 
along this direction to localize the one-magnon excita- 
tion along this direction in space. If the dispersion is 



completely flat, the magnon can be localized in a finite 
region of the lattice. This localized excitation can have 
N/n s different positions. Now it is clear that one can 
construct further local excitations. There will be no in- 
teraction with the other excitations as long as they are 
sufficiently well separated in space and therefore each ex- 
citation will have the same energy. In this manner, one 
obtains n-magnon excitations for n < n max whose energy 
is exactly n times the one-magnon energy. Due to the ab- 
sence of attractive interaction, it is plausible that these 
excitations are also the lowest n-magnon excitations. A 
proof of this statement for s = 1/2, arbitrary A > 
and all Jij equal will be given elsewhere |l2| - below we 
will report numerical evidence for the occurrence of this 
effect for several models. The essence of the above argu- 
mentation is that the ground state energies in the 1-, 2-, 
. . . , rt max -magnon spaces depend linearly on the number 
of magnons, i.e. on the total magnetic quantum number. 
Hence the total S z in the ground state goes directly from 
Ns — n max to the saturation value Ns when increasing 
the magnetic field. In terms of the magnetization curve 
m(h) = S z (h)/(Ns), this implies that there is a jump 
5m = n max / (Ns). If one band is completely flat, the sys- 
tem can support a macroscopic amount of independent 
magnons n max ~ N and one obtains a macroscopic jump 
just below saturation. 

To be more precise, denote the region of localization 
of the magnon state by L. Then the coefficients a; in (j^) 
are different from zero only for sites I G L. The local 
one-magnon state is completely decoupled from the rest 
of the lattice R and the eigenstate | ^f) can be written as a 
product \\V) — \^l)\^r) of a local part L and the rest R. 
\*&l) is the local magnon state and \^r) is the vacuum 
state. The coefficients <Zj vanish for i G R in the one- 
magnon state (|^), i.e. at =^ Vi G L and a, = Mi G R. 

The necessary and sufficient condition for decoupling 
of the local state from the rest R is 

^a ; J ;fe = V k£R. (3) 
leL 

The Hamiltonian (Q) can be divided into three parts 

H = H L + H L _ R + H R (4) 

with L being the part of the lattice where one magnon is 
localized and R the rest. The first term Hl is the local 
part of the Hamiltonian with — Ji t i 2 and G L. 
whereof is the lowest eigenstate. The second term 
Hl-r is the coupling of the local part to the rest of 
the lattice with J^j = Jik and / G L, k G R. Jik must 
satisfy condition (]3|) . The rest of the Hamiltonian which 
is not connected with the local part is Hji with = 
Jk ± k 2 an d ki,k2 G R. Hl-r together with condition 
creates the frustration, therefore it seems that the 
magnetization jump described here is restricted to highly 
frustrated lattices. The simplest realizations of such a 
Hamiltonian are rings connected only by triangles. 



2 




FIG. 1: Kagome spin array, which hosts at least N/9 indepen- 
dent magnons in a y/3 x y/S structure if it consists of N spins. 
The circles mark hexagons where the independent magnons 
can be localized. The structure of one localized magnon is 
indicated by the signs around the bold hexagon which corre- 
spond to coefficient ai = ±1 with which a spin-flip contributes 
at each site I. Together with the surrounding triangles the lo- 
cal state satisfies equation (^|). 

In the kagome lattice Jl], [?J as a typical example of 
flat one-magnon dispersion wo(k) = h— 2sJ(l + 2A) the 
magnon can be localized around a hexagon (see Fig. |l|). 
Choosing the coefficients a; = (—1)' with I numbering 
the sites around a hexagon ensures that is an exact 
eigenstate for the hexagon (as illustrated by the bold one 
in Fig. flh. The triangles around the hexagon fulfill con- 
dition (ph and therefore the magnon on the hexagon is 
decoupled from the rest of the lattice, showing that this 
state is also an eigenstate of the whole kagome lattice. 
Further magnons can be put on the lattice without dis- 
turbing existing excitations. This can be repeated until 
every third hexagon is excited as shown in Fig. |l|. As a 
consequence we have a macroscopic magnetization jump 
dm having its maximal value Sm = 2/9 for the extreme 
quantum case s = 1/2. 

A proof that these states are the lowest eigenstates in 
the corresponding sectors is given in [|l2| for s = 1/2, 
arbitrary A > and all equal. Here we present nu- 
merical evidence for this statement. In Fig. || exact di- 
agonalization results are shown for finite systems with 
periodic boundary conditions for the s = 1/2 Heisenberg 
antiferromagnet on the kagome lattice. Only lattices with 
N being a multiple of 9 (three unit cells) are presented, 
which fit to the m = 7/9 state corresponding to Fig. [|. 
The jump to saturation can be easily seen in this figure. 
We have also computed curves on smaller clusters. They 
agree with those presented in M, || . If the boundary con- 
ditions do not fit to the y/3 x V3 state or if the cluster 
is too small, the jump may show finite-size effects. 

We emphasize that the existence of this jump is quite 
independent of several details of the system. First, one 
can construct iV/9 independent magnons for any s, lead- 
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FIG. 2: Magnetization curves of an s — 1/2 Heisenberg an- 
tiferromagnet on a finite kagome lattice for N = 27, 36, 45, 
54 and A = 1. For N = 45 and iV = 54 the curve starts just 
below the magnetization jump. It can be clearly seen that 
the jump to full magnetization has no finite size dependence. 



cal magnon excitations. The lowest two out of the four 
magnon bands are indeed degenerate and completely flat: 
w (k) = wx{k) = h- s(J + J')(l + 3A). We expect a 
macroscopic jump of 5m > l/(12s) for all J, J' > 0. 
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ing to a jump of height 5m = l/(9s). Indeed, the results 
for the s — 1 kagome lattice presented in |?j show a jump 
of the expected height for the given cluster sizes. Second, 
introduction of an AJf Z-anisotropy A ^ 1 does not af- 
fect the crucial properties of the one-magnon dispersion 
and therefore one expects the magnitude of the degener- 
acy and the associated jump to be independent of A. In 
the exact diagonalization results for s = 1/2 || jumps 
of identical size are indeed observed for A = 0, 1 and 
2.5. Third, the argumentation remains also unchanged 
if one generalizes to different coupling constants in the 
triangles pointing up and down (see Fig. |l|) ||. The de- 
generacy is (partially) lifted only if coupling constants 
are changed such that they become different around one 
triangle. The jump therefore seems to be very stable not 
only in the kagome lattice but also in the other mod- 
els to be discussed next where similar arguments can be 
applied. 

Another 2D example for completely flat one-magnon 
dispersion wo(k) = h — 2s J(l + 3 A) is the checkerboard 
lattice |n|, a 2D variant of the pyrochlore lattice. In this 
case, localized magnon excitations live around a square 
without diagonal interactions, again with coefficients 
ai = (—1)'. The magnetization jump is 5m — l/(8s). 
We have verified the predicted degeneracy and associ- 
ated macroscopic jump numerically for the checkerboard 
lattice with s = 1/2 and A = 1. 

Completely flat bands can also be found in dimen- 
sions different from two. For example, the generalized 
pyrochlore lattice in three dimensions with two differ- 
ent coupling constants J, J' (see e.g. Q) gives rise 
to the high frustration necessary for a decoupling of lo- 



FIG. 3: One dimensional quantum antiferromagnets with a 
magnetization . j um p to the full moment m = 1: a) The saw- 
tooth chain [ jl5[ p"q ], b) the kagome like chain of and c) 
the kagome like chain of . The position of a local magnon 
is marked by thick lines with coefficients a; attached. 

Even in one dimension, one can find systems with a 
flat dispersion: Some examples are shown in Fig. || to- 
gether with the structure of the localized magnon exci- 
tations. For the generalized sawtooth chain [ll| of 
Fig. ||a), the lowest magnon branch is completely flat for 
J 2 = \/2(l + A) Ji and all s. Note that this example sat- 
isfies (y) with more complicated coefficients ai than en- 
countered previously. The lowest magnon branch for the 
one-dimensional kagome variant |l7j shown in Fig. |j^b) is 
also completely flat w (k) = h - 2sJ(l + 2A). Fig. |c) 
shows another variant of a kagome chain Here, the 
state indicated by the bold hexagon is an eigenstate for 
J 2 = (2A + l)Ji/(A + l) with w {k) = h-2sJ 1 (l + 2A). 
We have checked numerically that for s = 1/2 and A = 1, 
a jump of size 5m — 1/2, 1/3 or 1/5 exists in case a), 
b) or c), respectively. As an example, Fig. || shows the 
m(h) curve of the model of Fig. ||c) with s = 1/2, A = 0, 
and Ji — Jx (to ensure a flat dispersion). The jump of 
height 5m = 1/5 can be seen clearly (compare also the 
inset). Furthermore, one can see several plateaus in the 
magnetization curve. This suggests that the same condi- 
tions which give rise to the jump also favor the formation 
of magnetization plateaus. In particular, our numerical 
data always show a plateau preceding the jump. How- 
ever, this is beyond the scope of the present Letter and 
needs further investigations. 

So far, we have discussed cases with completely fiat 
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FIG. 4: Magnetization curves of the s = 1/2 YY-kagome 
chain (A = 0) of Fig. |c) with Ji = ,h = J for N = 20, 30, 
40, 50. Inset: Region just below the jump to saturation. 



dispersion. However, in more than one dimension it is 
also possible that the lowest magnon branch has a flat 
dispersion only in some, but not all directions. This 
has been noticed previously for the J\ — Ji model on 
the square lattice antiferromagnet with J 2 = Ji/2 [ jl9| . 
The generalized checkerboard lattice jl3| is another 2D 
model which has flat directions if the couplings J' along 
the diagonals are larger than those for the square lattice 
J. In these cases, magnons remain delocalized along one 
dimension and thus the localization region L is a one- 
dimensional object. For the J\ — Ji model, L extends 
along one of the axes of the square lattice while for the 
checkerboard model it extends along one of the diago- 
nals. In both cases, ai = (—1)' satisfies the condition (||) 
and yields localized magnon excitations. However, now 
the number of possible magnon excitations is bounded 
by the linear extent of the lattice and therefore their 
density tends to zero in the thermodynamic limit. Con- 
sequently, the jump is only observable in finite systems 
(this has been checked for both models with s = 1/2 and 
A = 1) and vanishes in the thermodynamic limit. Argu- 
ments similar to those used in [^| then indicate that all 
derivatives d l h{m)j dm % with i > 1 vanish at the point 
of full magnetization m = 1. This is in contrast to pfj 



lattices that for magnetic fields near saturation non- 
interacting magnons can condensate into a single-particle 
ground-state leading to a macroscopic jump in the mag- 
netization curve. This is a true quantum effect which 
vanishes if the spins become classical (s — > oo). 

Since this effect is generic in highly frustrated magnets, 
we are confident that a realization of some model dis- 
cussed in the present Letter can be found with sufficient 
small s and J to make the saturating field accessible. The 
dynamics of the lattice may be relevant in such a com- 
pound, but one will have to see whether this strengthens 
or weakens the anomaly predicted at the saturating field. 
Our analysis of deformations of the coupling constants 
at least indicates that the low-temperature magnetiza- 
tion curve will show an unusually steep rise even if the 
geometry deviates from the ideal structure. It is inter- 
esting to note that a similar jump also appears in certain 
molecular magnets fl2|| where the magnetic ions form an 
icosidodecahedron (like {Mx>72Fe3o}), a cuboctahedron or 
similar frustrated structures. 
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a square root was interred m m(h) tor the J\ — Ji model 
with Ji = Ji/2. Since the transition appears to become 
continuous for N — > 00, a non-analytic functional de- 
pendence in him) like him) ~ exp(— 1/(1 — to)) appears 
more plausible. In any case, also in such a situation the 
magnetization curve will be exceptionally steep close to 
saturation, a fact which should also be detectable in high 
field experiments on materials realizing such models. 
In summary, we have shown for certain frustrated spin 



[15. 
[16. 

[17; 

[18 

[19 
[20 



P. Lecheminant, B. Bernu, C. Lhuillier, L. Pierre, P. 
Sindzingre, Phys. Rev. B 56, 2521 (1997). 
Ch. Waldtmann, H.-U. Everts, B. Bernu, C. Lhuillier, P. 
Sindzingre, P. Lecheminant, L. Pierre, Eur. Phys. J. B 2, 
501 (1998). 

F. Mila, Phys. Rev. Lett. 81, 2356 (1998). 
H. Kageyama, K. Yoshimura, R. Stern, N.V. Mushnikov, 
K. Onizuka, M. Kato, K. Kosuge, CP. Slichter, T. Goto, 
Y. Ueda, Phys. Rev. Lett. 82, 3168 (1999). 
K. Onizuka, H. Kageyama, Y. Narumi, K. Kindo, Y. 
Ueda, T. Goto, J. Phys. Soc. Jpn. 69, 1016 (2000). 
T. Momoi, K. Totsuka, Phys. Rev. B 62, 15067 (2000). 
K. Hida, J. Phys. Soc. Jpn. 70, 3673 (2001). 
D.C. Cabra, M.D. Grynberg, P.C.W. Holdsworth, P. Pu- 
jol, Phys. Rev. B 65, 094418 (2002). 
M.E. Zhitomirsky, A. Honecker, O.A. Petrenko, Phys. 
Rev. Lett. 85, 3269 (2000). 

M.E. Zhitomirsky, Phys. Rev. Lett. 88, 057204 (2002). 
C. Gerhardt, K.-H. Mutter, H. Kroger, Phys. Rev. B 57, 
11504 (1998). 

J. Schnack, H.-J. Schmidt, J. Richter, J. Schulenburg, 
Eur. Phys. J. B 24, 475 (2001). 



B. Canals, P hys. Rev. B (to be published), cond 



mat/0102233 



B. Canals, C. Lacroix, Phys. Rev. Lett. 80, 2933 (1998). 
T. Nakamura, K. Kubo, Phys. Rev. B 53, 6393 (1996). 
D. Sen, B.S. Shastry, R.E. Walstedt, R. Cava, Phys. Rev. 
B 53, 6401 (1996). 

Ch. Waldtmann, H. Kreutzmann, U. Schollwock, K. 

Maisinger, H.-U. Everts, Phys. Rev. B 62, 9472 (2000). 

P. Azaria, C. Hooley, P. Lecheminant, C. Lhuillier, A.M. 

Tsvelik, Phys. Rev. Lett. 81, 1694 (1998). 

A. Honecker, Can. J. Phys. 79, 1557 (2001). 

M.S. Yang, K.H. Mutter, Z. Phys. B 104, 117 (1997). 



